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Abstract 

We evaluate the Casimir vacuum energy at finite temperature associated with 
the Maxwell field confined by a perfectly conducting rectangular cavity and show 
that an extended version of the temperature inversion symmetry is present in this 
system. 
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1 Introduction. 



Temperature inversion symmetry occurs in the free energy associated with the Casimir 
effect [m at finite temperature and it is hnked to the nature of the boundary conditions 
imposed on the quantum vacuum oscillations. As shown by Ravndal and Tollefsen 
temperature inversion symmetry holds for massless bosonic fields and symmetric boundary 
conditions and also for massless fermionic fields and antisymmetric boundary conditions. 
Temperature inversion symmetry was initially studied by Brown and Maclay who wrote 
the scaled free energy associated with the Casimir effect at finite temperature as a sum of 
three contributions, namely: a zero temperature contribution, that is, the Casimir energy 
density at zero temperature, a Stefan-Boltzmann contribution proportional the fourth 
power of the scaled temperature ^ := Td/n, and a non-trivial contribution. The three 
contributions to the scaled free energy can be combined into a single double sum: 

m) = - J- f f (^"-a' (1) 

where the term corresponding to n = m = is excluded. Setting n = and summing 
over m with m 7^ 0, we obtain the Stefan-Boltzmann limit — 7r^^''/45. On the other hand, 
setting m = and summing over n with n 7^ 0, we obtain the zero temperature Casimir 
term — 7r^/720. This function has the following property: 

(20 V (1/40 = / (0 , (2) 

which is the mathematical statement of the temperature inversion symmetry. It was also 
shown by Gundersen and Ravndal that the scaled free energy associated with massless 
fermions fields at finite temperature submitted to MIT boundary conditions satisfy the 
relation given by Eq.(|^) and therefore exhibts temperature inversion symmetry. Tadaki 
and Takagi [Q have calculated Casimir free energies for a massless scalar field obeying 
Dirichlet or Neumann boundary conditions on both plates and found this symmetry. For 
the parallel plate geometry and mixed boundary conditions it is possible to circunvent 
the restrictions found by and discuss temperature inversion symmetry by reconizing 
that with respect to the evaluation of the free energy this arrangement is equivalent to 
the difference between two Dirichlet (or Neumann) plates [^]. In the case of a massless 
scalar field at finite temperature and periodic boundary conditions, it is possible to show 
that the partition function, and consequently the free energy, can be written in a closed 
form such that the temperature inversion symmetry becomes explicit [0. Here we shall 
show that for the case of the Maxwell field confined by a perfectly conduting rectangular 
cavity this symmetry is also present, provided that we generalize the transformations that 
relate the high and the low temperature regimes. 
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Throughout this letter we employ units such that Boltzmann constant ks, the speed 
of light c and h = h/2TT are set equal to the unity. 



2 Evaluation of Helmholtz free energy 

Helmholtz free energy for the confined Maxwell field within a perfectly conducting rect- 
angular box can be evaluated by means of the generalized zeta function technique [0], 



IG]. First we introduce the global zeta function which is defined by 
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where f3 = 1/T , the reciprocal of the temperature T, is the periodic length along the Eu- 
clidean time direction and a mass scale parameter /x was introduced in order to keep Eq.(|^) 
dimensionless; u^^^y are the eigenvalues associated with the Euclidean time-independent 
modal equation 

- Ax V{n} (x) = <^{n} (x) , (4) 

where x = (x, y, z) and Ax is the Laplacian operator. The Helmholtz free energy function 
F {(3) for bosonic fields can be obtained from the knowledge of the global zeta function 
through the relation 

1 a r 
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The eigenfrequencies associated with the allowed electromagnetic normal modes within a 
perfectly conducting rectangular box of linear dimensions a, 6, c are given by 
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where l,m,n G {1,2, . . .}. For l,n,m ^ there are two possible polarization states and 
for / = 0orm = 0orn = there is one polarization state only. Eigenfrequencies for which 
three or two of the integers l,m,n are simultaneously zero are not allowed. It follows that 
the generalized zeta function for the case in question reads 
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If we now separate the terms corresponding to the zero temperature sector by setting 
p = 0, then it will be easily seen that we can rewrite Eq.(0) formally as a sum of Epstein 
functions (see the Appendix) which reads 



TT 



C (s) = {s; l/a\ l/b\ l/c^ f) + 23^3 (s; l/a^, l/fo^ i/c^ 

+23^3 [s- l/a\ l/h\A/l3^) + 22^2 (s- l/a\ 1/6^ 
+2=^^3 {s- l/a\ l/c\ A/ (3^) + 22^2 (s; l/a\ l/C 
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In order to regularize Eq.(D it is convenient to rewrite it as a sum of generalized zeta 
functions and this can be accomplished by making use of Eqs. (^) and ( [26| ) in the 
appendix. Equation (|]) then takes the form 
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Finally, we can regularize Eq.(^ by applying the reflection formula, Eq.(^) on all terms 
except the last one which is already regularized thus obtaining 
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where Cj? (2-2) is the Riemann zeta function and we have also made use of the fact that 
E\ (2; a) = a~^C/j(22;). Notice that all terms on the R.H.S. of Eq. (p!0|) are zero when 



s — i> except the last one, that is, C (0) 7^ since (r{0) = — (1/2) ln(27r). This will 
give rise to a scale-dependent term absent in previous calculations related to this problem 
1^. By making use of Eq.( |^) we see that the Helmholtz free energy for the confined 



electromagnetic field confined within a perfectly conducting rectangular box at finite 
temperature is 
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The presence of a scale-dependent term is a leftover from the renormalization procedure 
carried on via the generalized zeta method and introduces an ambiguity in the Helmholtz 
free energy related to the problem in question. This ambiguity is introduced whenever C{s) 
is not zero at s = 0, see for instance[jlO|- It was shown by Myers 0] that this ambiguty 



is a legitimate feature of the generalized zeta function regularization method, whenever 
the theory is massive or/and interacting and/or some dimensions are compactified. As 
consequence the usual mode sum method and the generalized zeta function method are 
not always equivalent. However, the usual mode sum method can be modified so as to 
include the complexities of a scale dependence |12[. At zero temperature and for the 



simple geometrical cavity that we are considering here scale- dependent terms are not 
expected. This is confirmed by the independent calculations due to Lukosz , Ruggiero 
et al [^, Ambjorn and Wolfram |jl5[ and the present authors. Here it is clear that the 
scale dependence is introduced by the generalized zeta function regularization method due 
to the imposition of periodic conditions on the Euclidean time direction. For our purposes 
the ambiguity in the free energy can be solved by the additional physical requirement that 
in unconstrained space, that is, for a very large rectangular box, the only surviving term 
must be the Stefan-Boltzmann term. It is seen then that /i should be equal to \/2t: / (3. 



3 Temperature inversion symmetry 

The explicit verification is of temperature inversion symmetry in the case of the system in 
question is complicated by the fact that in addition to the inverse temperature parameter 
(3 we have to deal with three other length parameters, namely, a, h and c, the measures 
of the sides of the rectangular cavity. At finite temperature and periodic or antiperiodic 
conditions along one spatial direction, or Dirichlet, or Neumann plane surfaces located 
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perpendicularly to one spatial direction only two characteristic lengths are involved, a 
feature that rends this verification easier. Neverheless, by making recourse to a simple 
trick we shall show that temperature inversion symmetry is also present here. 

First notice that the denominator of the first term on the R.H.S of Eg. (pi]) can be 
rewritten as a^/^ + b'^m? + c^n^ = {a?q^ + 6^r^ + c^t^), where {g, r, t] is a sequence of 
three integers with no common factor and k is the common factor of {l,m.n}. For the 
sequence {g, r, t} we define a characteristic length d{q^r,t} by 

d\^^^^^ := + ftV^ + cH^. (12) 

Let us also define the dimensionless variable 

4{q,r,t} •- ^ • [i-'i) 

The remaining terms on the R.H.S. of Eq.(^) can be treated in a simpler way. For 
example, for the second term we define 

^. (14) 

Then making the replacement: Y.i,m,n=-oo ~^ Y.{q,r,t}Y.t=-oo^ f^^^ energy density 
corresponding to Eq.(^) can be written as 
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The zero temperature limit of Eq.(|T^) is obtained by reconizing that when /5 — oo 
the dimensionless parameters ^{q,r,t},^a,^b,^c — ^ and the surviving terms in Eq. (p!5|) 
correspond to j9 = 0, as the reader can easily verify. Going back to our initial set of 
indexes l,m,n we obtain 

^ a6c +2? 1 vr /I 1 1\ 
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which is in perfect agreement with the result obtained by Lukosz [O] and also by Ruggiero 



et al. 14 



In order to introduce temperature inversion symmetry for this problem let us now 
define the dimensionless functions 



167r2 

K,p= — oo 



and 

1 oo a2 



with similar definitions for JF^ (^fe) and JF^ (^c)- Therefore, in terms of these functions we 
can write the free energy density as 

V ih ^U} ^ aV ^ bV ^ cV ^ ^''^ 

where V = abc. It is not hard to see that the functions ^{q,r,t} {^{q,r,t}j exhibit the 
following property 

^{q,r,t}^{q,r,t} ^ ^ = -^{g,r,t} {^{q,r,t}) ■ (20) 

In the same way 

C fa = fa i^a) , (21) 

Equations (^) and (pTD plus two similar equations for the functions J-'h {^b) and J-'c (^c) 
describe temperature inversion symmetry for the case in question, that is, all terms in 
Eq.(|19D can be inverted by making use of these formulae. 

In the very high temperature limit we expect the leading contribution to be the Stefan- 
Boltzmann term, vr^/ (45/3'^). Now we argue that our transformations are applied in order 
to generate from the unique Stefan-Boltzmann density term, an infinite number of terms 
which must be added at zero temperature, and inversely, each term at zero temperature 
goes to the unique Stefan-Boltzmann term. In this form, that is, as far as the energy 
density is concerned, we can apply the temperature inversion symmetry transformations to 
several Casimir systems, including of course the ones previously analyzed in the literature. 
The specific form of the transformations depends on the particular system at hand, but 
the procedure is the same. 
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Now, we can easily check that J^{q,r,t} (0) = — 7r^/720, and that J-'a (^a) = i^b) 
(^c) = 7r/48. Then, making use of Eqs. pO|) and (pT)) we write 



i{q,r,t)^{q.rfy (oo) = - — , (22) 



720' 

and 

fool = Tu fool = ^2 fool = 

48 



il (oo) = (oo) = -Fc (oo) = ^. (23) 



Taking these resuhs into Eq.(|T9D we obtain 



F (a, 6, c, /5 ^ 0) ^ + __ (a + 6 + c) , (24) 



which is in agreement with |T5[ with respect to the leading terms in the high temperature 
approximation. As expected, the Stefan-Boltzmann term is the leading term in this limit. 

It is convenient to rewrite Eq.(|^), or its equivalent Eq.(|T^), with the zero tempera- 
ture terms and the Stefan-Boltzmann term separated from the non-trivial temperature 
corrections, 

F TT^ vr / 1 1 1 
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The last two lines of Eq. (|25D represent the non-trivial temperature-dependent corrections. 
The first term in the third line is equivalent to a set of conducting plates with each pair 
of plates characterized by a distance d^q^r,t}, which means that we can take advantage of 
results already established in the literature regarding high and low temperature approxi- 
mations for a pair of conducting plates. 

4 Conclusions 

In this letter we have employed generalized zeta function techniques in its global version 
to derive the finite temperature vacuum energy associated with an electromagnetic field 
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confined within a perfectly conducting rectangular box. We have shown that Helmholtz 
free energy for the problem in question is scale dependent, a feature that was not present 
in previous calculations concerning this problem at zero temperature [0], ||14[ and finite 
temperature [|1^. The scale dependence is a remainder that tell us that the theory was 



renormalized by the way of the generalized zeta function method. We have also shown 
that the concept of temperature inversion symmetry can be extended so as to include 
the elctromagnetic field confined by perfectly conducting rectangular cavity. The main 
point here is the recognition that the concept of temperature inversion symmetry can be 
extended in the sense that it must be applied separately to the different terms comprising 
the non-trivial part of the free energy density, with each term behaving differently under 
this symmetry. This is the content of Eqs. (^OD and (pT]) . The most important feature of 
temperature inversion symmetry is that it allows us to establish a relationship between the 
zero and low-temperature sector of the Helmholtz free energy and the high-temperature 
one. This can be very useful since in general it is easier to obtain low-temperature 
expansions. In particular, it is possible to relate the zero temperature Casimir effect to 
the Stefan-Boltzmann term in a straightforward way. 



Appendix 



Multidimensional homogeneous Epstein function (see also Elizalde et a/.|10| for a mod- 
ern presentation) are defined by 



Eat (z; oi, ...Oat) = (ainl + a2nl-\ hOArn^) , (A.l) 

for > and oi, 02, ...Oat > 0. Generalized zeta functions. Aat (-2; oi, 02, •••Otv) are 
defined by §,[[101 

+00 / _^ 
An {z;ai,a2, ...aN) = E (ainl + a2nl -\ hOArn^) , (A.2) 

ni,n2,..npf=—oo 

for > N/2 and ai,a2, ...aN > 0. The prime means that the term corresponding 
to n-i = 77-2 = ... = ?T-Af = must be excluded from the summation. A useful reflection 



formula reads envolving these functions is |T0 
An {z; Oi, 02, ...ajy) 



TT 2 
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Relationships between generalized zeta functions and Epstein functions can be derived 
from the very definition given by (^). For example, 

A2 {z\ ai, aa) = 1^E2 {z\ ai, aa) + {z\ a^) + {z\ 02) {AA) 
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and 



A4 {z; ai, a2, as, a^) = 2'^i?4 {z; Oi, 02, as, 04) + 2^i?s (-2; ai, a2, as) 

+2^£'3 (^; ai, a2, 04) + 2^£;s {z; ai, as, a4) + 2^£'2 {z; ai, a2) + 2^£'2 {z; ai, as) 

+2^£;3 {z; 02, as, 04) + 2^£;2 {z; 02, as) + 2^£;2 {z; ai, a4) + 2'^E2 {z; aa, a4) 

+2'^E2 (z; as, a4) + 2Ei {z; ai) + 2E1 {z; aa) + 2£;i {z; as) + 2£;i (z; a4) , (A.B) 

of which use was made in this letter. 
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